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The level curvature distribution function is studied beyond 
the random matrix theory for the case of T-breaking pertur- 
bations over the orthogonal ensemble. The leading correction 
to the shape of the level curvature distribution is calculated 
using the nonlinear sigma-model. The sign of the correction 
depends on the presence or absence of the global gauge invari- 
ance and is different for perturbations caused by the constant 
vector-potential and by the random magnetic field. Scaling 
arguments are discussed that indicate on the qualitative dif- 
ference in the level statistics in the dirty metal phase for space 
dimensionalities d < 4 and d > 4. 
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Since the seminal work by Thouless [1] the sensi- 
tivity of spectrum {En} of disordered conductors to a 
small twist of the boundary conditions 5'(p, z = 0) = 
e"''4'(p, z = L) is considered as a powerful tool to probe 
the space structure of eigenfunctions and distinguish be- 
tween the extended and the localized states. More pre- 
cisely, the quantity Kn which is now referred to as the 
"level curvature", was introduced in Ref.[l] in order to 
describe this sensitivity quantitatively: 
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where the mean level spacing A — [vL"^)^^ is related 
with the mean density of states v and the size of the 
rf-dimensional sample L. 

In complex or disordered quantum systems, the quan- 
tity Kn fluctuates over the ensemble of energy levels {En} 
or, for a given level, over the ensemble of realizations of 
disorder. The typical width of the distribution of level 
curvatures P{K) is of the order of the dimensionless con- 
ductance g = Z3/(L^A), where D is the diffusion coeffi- 
cient [1] . Thus studying this distribution one can obtain 
an information about the transition from metal to insu- 
lator with increasing the disorder (Anderson transition). 

In the recent works [2] it has been shown that the 
distribution of the fluctuating quantity Kn is a partic- 
ular example of parametric level statistics, i.e. statistics 
of spectral responses of the system to any perturbations 
proportional to some parameter if. 

A remarkable property of the parametric level statis- 
tics [2] similar to that of the usual level statistics [3,4], 
is that in a certain limit they are universal for all clas- 



sically chaotic and disordered systems and can be de- 
scribed by the random matrix theory (RMT) of Wigner 
and Dyson [3,4]. For disordered systems considered here 
this limit coincides [5] with g ^ oo. For chaotic sys- 
tems the same role is played [6] by the ratio g — 71 /A, 
where 71 is the first non-zero mode in the spectrum of 
the Perron-Frobenius operator that describes the chaotic 
behavior of the corresponding classical system. In partic- 
ular, for the time-reversal-invariant system without spin- 
dependent interactions (orthogonal ensemble) the distri- 
bution of level curvatures, Eq.(|l|), was found [7,8] in this 
limit to have the form: 
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where {\K\)rmt — 2g is the average modulus of the 
level curvature. Further study [9] showed that the form, 
Eq.(H), remains the same if weak localization is taken 
into account, only the dimensionless conductance in the 
expression for {\K\) is changed appropriately. 

The form, Eq.(||), is universal. It does not depend, 
e.g. on details of the system and the perturbation. The 
only important point is that the system is T-invariant 
with = 1 and the perturbation breaks this invari- 
ance. However, the T-breaking perturbations can be of 
two different types. One of them (we call such perturba- 
tions " global" ) is represented by the twist of the bound- 
ary conditions or by the magnetic flux in the problem of 
persistent current [10]. In this case wave functions should 
have a significant amplitude on both boundaries z — 
and z ~ L in order the energy spectrum En be sensitive 
to the phase (p. Apparently, it is not so for the localized 
states, for which the distribution function P(k) shrinks 
dramatically and is approximately of the logarithmically- 
normal form [11-13]. 

Another type of T-breaking perturbations is repre- 
sented by the magnetic impurities or random magnetic 
fluxes. We call them "local perturbations", since they 
probe wave functions locally. Even for a localized state 
there should be a significant spectral response to chang- 
ing the strength ip of such a local perturbation. So, one 
would expect neither any dramatic contraction nor the 
logarithmically-normal form of the distribution function 
P{K) in this case. 

Thus we conclude that there is no universality of the 
distribution function P{k) in the Anderson insulator. It 
is reasonable to assume that the universality is broken 
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already at small but finite values of the parameter \/g 
when the wave functions are still extended but have a 
pronounced structure. In the presence of such a struc- 
ture, the system cannot be described any more by the 
basis-invariant ensemble of random matrices which is the 
main assumption of the random-matrix theory of Wigner 
and Dyson [3,4]. In order to use the curvature distribu- 
tion function as a tool to probe the space structure of 
wave functions one should go beyond the random matrix 
theory. 

In the present Letter we address the problem of cor- 
rections to the shape of the distribution function P{k). 
A general approach to calculate such corrections using 
the nonlinear supersymmetric sigma-model [5] has been 
suggested in [14] and applied to distributions of differ- 
ent quantities [9,14,15]. It is based on the perturbative 
consideration of the non-zero diffusion modes which are 
integrated out to produce corrections to the zero-mode 
supersymmetric sigma-model [5] . The latter is then inte- 
grated exactly using the Efetov's parametrization [5]. 

Before going into details of the calculations we would 
like to formulate the main result and its consequences. 
The main result is the correction to the level curvature 
distribution function P{k) for the case of T-breaking per- 
turbations over the orthogonal ensemble: 
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Here q = {qi, ...Qd}, where qi = 27rnj, (n^ = 0, ±1, ±2...) 
in the case of the periodic boundary conditions (for an 
unperturbed system) considered in this Letter. 

A remarkable feature of Eqs.(^)-(^) for < 4 is that 
the sign of the correction is different for global (case I) 
and local (case II) T-breaking perturbations. This re- 
flects the principal difference between them which we 
discussed above. The positive sign of the correction at 
small K for the case of the global perturbations has been 
observed in the recent numerical simulations [12] on the 
3D Anderson model. It implies the tendency towards 
the formation of a sharp peak in P{K) at if = which 
should be present in the insulator for the case of global 
perturbations and absent in the case of local ones. 

Next we note that the sum in Eq. (^ is convergent only 
for d < A. For d > 4 there is a problem of the upper mo- 
mentum cut-off that cannot be done accurately within 
the framework of the nonlinear sigma-model approach. 
It turns out that for global perturbations the coefRcient 
C{d) changes sign precisely at d = 4 where the difficulty 
with the divergent sum is first encountered. This coinci- 
dence could hardly be accidental. It may indicate on the 
qualitative difference in the structure of wave functions 
in disordered systems for d < A and d > 4, especially near 



the Anderson transition where the parameter 1/g^ is not 
very small. Meanwhile the problem of the Anderson tran- 
sition in higher dimensions is becoming a physical issue 
because of some models [16] that map the problem of 
disordered systems of interacted electrons onto the An- 
derson model in higher dimensions. 

Now we outline the derivation of the announced re- 
sults. The first step is to express the level curvature 
distribution function P{K) in terms of the two-level cor- 
relation function R{uJ, ip). This is done in the same way 
as in Ref.[17]: 
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The second step is to represent the two-level correla- 
tion function in Eq.(|^) in terms of the supersymmetric 
nonlinear sigma-model which has been derived in Ref. [5] 
from the model of non-interacting electrons in a random 
potential. Using the results of Ref. [5] we arrive at: 
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In this equation Q{r) is the 8x8 supermatrix field of 
a certain symmetry which obeys the constraint — 
1, and F^[Q; K,j±] = F[Q; K] + F[Q; j+] + F[Q; j_], 
where: 



f dr 

F[Q; j±] = j± J Str[P±PBQ] 



(7) 



Here S'tr[^] = ir[AFi^-] — ir[^Bs] stands for the super- 
trace, P± are the projectors onto the retarded (RR) or 
advanced (AA) sectors and Pb is the projector onto the 
boson-boson sector of the superfield Q (see Ref. [5] for de- 
tails and definitions). An equivalent representation has 
been used in Ref. [9] . 

The main information about the underlying physics is 
contained in the action F[Q] K]. For global perturba- 
tions it takes the form: 

F^ [Q; if] = ^gL^ j Str (vQ -'-^[Q,t]\ + (8) 



where — (0, (f) is supposed to be directed along the 
z-axis; A = P+ — P_ , and r = ct^ ® P+ is proportional 
to the Pauli matrix az that breaks the symmetry in the 
quaternionic subspace of the superfield Q and reflects the 
T-symmetry breaking. 

For local T-breaking perturbations the linear in Lp term 
in Eq.(||) is absent. 

The representation, Eqs.(|6|)-(||), in terms of the field Q(r) 
contains all the spatial diffusion modes 7g = {D/L'^)c^. 
However, in doing the limit — > in Eq.(^ the main 
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role is played by the zero mode which corresponds to 
q = 0. At (/? = this mode does not cost any energy no 
matter how large are the components of the field Q in the 
non-compact boson-boson sector [5]. It is the arbitrary 
large amplitudes of the zero mode components of the 
field Q that compensate the infinitesimal parameter ip in 
Eqs.(|)-(||) and lead to a finite resuh for P{K). Thus 
the space independent zero mode Qo niust be considered 
non-perturbatively. In the limit 5 ^ cxd all the non-zero 
modes can be neglected [5], and one arrives [9] at the 
RMT resuh, Eq.(|). For finite I/5 the non-zero modes 
should be also taken into account. However, all the non- 
zero modes can be treated perturbatively for 5 3> 1 to 
lead to some corrections to the zero-mode action. In 
order to obtain these corrections we have to separate zero 
modes from all other modes and then integrate over all 
the non-zero modes using a certain perturbative scheme. 
This is done by means of the transformation [14]: 



(9) 



The field Q{v) that must also obey the constraint = 1, 
is parametrized as follows [14]: 



Q = (1 - Wl2)k{l - W/2)- 
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where the supermatrix W^(r) does not contain the zero 
spatial mode. 

The transformation from Q{y) to new variables 
Qo,Wq^o involves [15] a non-trivial Jacobian J = 
^-Fj[w\^ To the first non- vanishing order in W it is given 
by [18]: 



(11) 



In order to integrate over W perturbatively, one has to 
expand the action, Eqs.(|7|),(p|),(pl]) in powers of W sep- 
arating the Vo-independent quadratic in W part Fo[iy] 
and consider all the rest M^-dependent terms as pertur- 
bations using the Wick's theorem. 

As a result of such calculations we have obtained the 
correction up to the order to the zero- mode action 
F[(5o; K,j±]. It is too lengthy to present it here and 
contains all possible linear and bi-linear combinations of 
vertices encountered in -F[Qo; K^j±] plus 4 new vertices. 

A remarkable feature of the Efetov's parametrization 
[5] that has not been mentioned so far is that in the 
limit of large "non-compact angles" A1.2 — > 00, the new 
vertices reduce to the product of the "old" ones. Thus 
up to the leading order of X\ 2 we found: 

Str[{KQof] - - [Str[kQof , (12) 
Str[{TQ^f] = -]^[Str[{TQ^f]]^ 

5ir[(TQo)'AOo] - -]^Str[{TQof] Str[KQol 

Str[{TQofP±PBQo] - -lstr[{TQo)^] Str[P±PBQo]- 



This circumstance simplifies the solution dramatically, 
since it allows to obtain all the corrections to P{K) just 
by differentiating of the RMT result, Eq.(^), over g and 
K. Since the K and ^-dependences of the RMT distribu- 
tion function Pwd{K) — {2g)^^PwDiK/2g) are closely 
related, the result can be formulated in terms of the g- 
derivatives alone. It is convenient to introduce the pa- 
rameter a = lug. Then the corrections to the distribu- 
tion function P{K) can be represented in the form: 
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The result, Eq.(|13[), holds for the global perturbations. 
For such perturbations VQ enters Eq. (||) as a " covariant 
derivative" VQ — ^[Q^t]. As a result, the linear in (p 
"cross" -vertex (p Str[QV Qt] appears in the action. It is 
this vertex that leads to terms in Eq. (|l3|) proportional 
to 1/d or lld{d + 2). For local perturbations such terms 
are absent, since the linear in ip vertex does not appear 
at aU [5]. 

The term with the second derivative in Eq. (p^) con- 
tributes to the correction to the shape of the distribution 
function P{k) and leads to the announced result, Eq.(||). 

The terms proportional to the first derivative can be 
absorbed into the RMT distribution function PwD{k) by 
changing the average (l^^l). The l/g correction to this 
value first found in Ref. [9] is nothing but the usual weak- 
localization correction to the quantity {\K\)b,mt — 2g. 
Thus up to the first order in the inverse "bare" conduc- 
tance 1/(7, the quantity ^([i^l) behaves exactly like the 
true dimensionless conductance g. 

Moreover, it turns out that all the 1/g^ corrections 
to P{K) proportional to (^ ^)^ cancel each other in 
the same way as in the dimensionless conductance g [20] . 
This is in a full agreement with the one-parameter scaling 
[21] extended to describe the distributions of mesoscopic 
fluctuations [19]. In the region where the one-parameter 
scaling works, the distribution functions of all quantities 
(measured in proper units) should depend only on one 
additional parameter - the average dimensionless con- 
ductance g that absorbs all microscopic parameters of 
the system and includes all quantum corrections propor- 
tional to the powers of ^ ^ . For the case of the level 
curvature distribution it means that P{K) = P(K;g), 
where g should obey the scaling equation [21] that does 
not allow [20] the term proportional to g~^{J2 ^)^ 
the perturbative expansion of g/g in powers of l/g. 

^From this viewpoint the corrections to P{K) propor- 
tional to g~^ ^ should be considered as the main 
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terms in the expansion of P{K; g) in powers of the true in- 
verse dimensionless conductance 1/g renormahzed by the 
locahzation effects rather than the second-order terms of 
an expansion in powers of the "bare" value of With 
due regard of these corrections the shape of the distri- 
bution function P{k) deviates from the RMT result and 
the average {\K\) is no more proportional to the average 
dimensionless conductance: 



\K\)=2g 
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where A{d) = 9/2 for the local perturbations and A{d) = 
(9/2 - 16/d + 36/d{d + 2)) for the global perturbations. 

In conclusion we would like to discuss a consequence 
of the form of corrections to different level statistics for 
the applicability of the RMT results. The typical relative 
correction [see e.g. [14,22]] has the same form as Eq.(p^) 
and is proportional to ^ (q^' suppose that 

the system is close enough to the Anderson transition 
(from the metal side) but it is still not in the critical 
region, so that irg ^ 1 but it is not very large. Then 
for c? < 4, where the sum in Eq.(p^) is convergent, the 
correction is still small, and the RMT can be considered 
as a good zero approximation. However, it is not the 
case for d > 4 if ng < (i/ro)''/^"^, where r^^ is the 
upper momentum cut-off. In the scaling picture of the 
Anderson transition [21] the dimensionless conductance g 
can be made arbitrary close to its critical value 1 by 
the fine tuning of the elastic scattering mean free path 
Z « at a constant L ^ Z ~ tq. Thus we conclude 
that the above condition is always possible to achieve for 
d > 4 at some strength of disorder Wd {L) which is less 
than the critical one Wc- For W > Wci the RMT results 
do not work even as a zero approximation for the level 
statistics. 

In summary, we have calculated the leading corrections 
to the shape of the level curvature distribution function 
P{K) beyond the random matrix theory. For T-breaking 
perturbations over the orthogonal ensemble, the correc- 
tions are not universal and depend qualitatively on the 
presence or absence of the global gauge invariance. In 
both cases the average |(|^|) shows deviations from the 
average quantum dimensionless conductance g. The cor- 
rections are in agreement with the one-parameter scaling. 
The consequence of the scaling hypothesis and the form 
of the corrections is a qualitative difference in level statis- 
tics for space dimensionalities above and below d = 4. 
Below d = A the RMT is a good zero approximation for 
level statistics in the metal phase, while for d > 4 there 
should exist a critical value of disorder Wd < Wc above 
which RMT fails to describe spectral correlations. 
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